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THE AVERAGE FIELD APPROXIMATION FOR ALMOST BOSONIC 

EXTENDED ANYONS 

DOUGLAS LUNDHOLM AND NICOLAS ROUGERIE 


Abstract. Anyons are 2D or ID quantum particles with intermediate statistics, inter¬ 
polating between bosons and fermions. We study the ground state of a large number N 
of 2D anyons, in a scaling limit where the statistics parameter a is proportional to N~^ 
when N ^ cx3. This means that the statistics is seen as a “perturbation from the bosonic 
end”. We model this situation in the magnetic gauge picture by bosons interacting through 
long-range magnetic potentials. We assume that these effective statistical gauge potentials 
are generated by magnetic charges carried by each particle, smeared over discs of radius 
R (extended anyons). Our method allows to take i? —>■ 0 not too fast at the same time 
as N ^ oo. In this limit we rigorously justify the so-called “average field approxima¬ 
tion”: the particles behave like independent, identically distributed bosons interacting via 
a self-consistent magnetic field. 
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1. Introduction 

In lower dimensions there are possibilities for quantum statistics different from bosons 
and fermions, so called intermediate or fractional statistics. Due to the prospect that such 
particles, termed anyons (as in anything in between bosons and fermions), could arise as 
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effective quasiparticles in many-body quantum systems confined to lower dimensions, there 
has been a great interest over the last three decades in figuring out the behavior of such 
statistics (see da dsi Eni Ea [391 iia [53] for extensive reviews). In one dimension, one can 
view the Lieb-Liniger model [55| as providing an example of effective interpolating statistics. 
Although initially regarded as purely hypothetical while at the same time offering substan¬ 
tial analytical insight thanks to its exact solvability, this system has now been realized 
concretely in the laboratory ED 03). Much less is known concerning fractional statistics in 
the two-dimensional setting, conjectured [T| to be relevant for the fractional quantum Hall 
effect (see dSlESI for review), and it is indeed a very challenging theoretical question to 
figure out even the ground state properties of an ideal 2D many-anyon gas, parameterized 
by a single statistics phase or a periodic real parameter a (with a = 0 corresponding 
to bosons and a = 1 to fermions). On the rigorous analytical side, some recent progress in 
this direction has been achieved in [341 [35l [36] where a better understanding of the ground 
state energy was obtained in the case that a is an odd numerator fraction. Numerous ap¬ 
proximative descriptions have also been proposed over the years, such as e.g. in [6| where 
the problem was approached from both the bosonic and the fermionic ends, with a harmonic 
trapping potential. Here, equipped with new methods in many-body spectral theory, we 
will re-visit this question from the perspective of a perturbation around bosons, i.e. in a 
regime where a is small. 

1.1. The model. One may formally think of 2D anyons with statistics parameter a as 
being described by an N-body wave function of the form 

- - TT .. ,XAr), 0,^ - arg^^— 

fA 

where ^ is a bosonic wave function, i.e. symmetric under particle exchange. This way 
T(xi, ...,Xj,...,Xk,.. .,xn) = ... ,Xk, ■ ■ ■ ,Xj,.. .,xn) 

and the behavior under particle exchange interpolates betweens bosons (a = 0) and fermions 
(a = 1). Of course the wave function is then in general not single-valued, and this descrip¬ 
tion is not easy to use in practice. To describe free anyons, one way out is to realize that 
acting on T with the free Hamiltonian 

N 

Y, (-A + v{x,)) 

i=i 

is equivalent to acting on the bosonic wave function T with an effective, a-dependent 
Hamiltonian. 

The Hamiltonian obtained in this way for N identical and ideal 2D anyons in a trapping 
potential H M+ reads (see [T5l [TH] [20l [23l [39l 021 [53] for review and [331 ESI [^ for 

recent mathematical studies) 

N 

Hn := ((P, + aAj? + V{x,)) (1.1) 

i=i 


where 


Pj = -iVj 
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is the usual momentum operator for particle j and (denoting {x,y)-^ := (—y, x)) 


E 


(Xj - Xk)-^ 
I , _ /y» , I 2 


( 1 . 2 ) 


is the (normalizecfl) statistical gauge vector potential felt by particle j due to the influence 
of all other particles. The statistics parameter is denoted by a, corresponding to a statis¬ 
tical phase 6*“^^ under a continuous simple interchange of two particles. In this so-called 
“magnetic gauge picture”, 2D anyons are thus described as bosons, each of them carrying 
an Aharonov-Bohm magnetic flux of strength a. 

We shall in this work assume 


a = 


/3 

1 


—)■ 0 when N —oo. 


(1.3) 


where /3 is a given, fixed constant. We consider (HH) as (formally) acting on the bosonic N- 
particle space Lgyj^(M^^), which together with the condition (11.31) means that we consider 
almost bosonic anyons. Note that if we simply took a —)• 0 at fixed A^, we would recover 
ordinary bosons at leading order. One could then only see the effect of the non-trivial 
statistics in a perturbative expansion, a route followed e.g. in [IHl O 091 SU [10] . However, 
if W —)• oo with fixed /3 as above, the anyon statistics has a leading order effect, manifest 
through a particular mean-field model with a self-consistent magnetic field of strength ~ /3, 
studied e.g. in [6l[50l[5Tl[l9l[52] and often called the average field approximation (see [^IT^ 
for review). Our aim in this work is to justify this description rigorously. 


Actually, the Hamiltonian dni) is too singular to be considered as acting on a pure 
tensor product however regular the function u of the one-particle space 

S) C L^(M^). We refer to [361 Section 2.1] for a discussion of the domain of which 

requires the removal of the two-particle diagonals from the configuration space One 

way to circumvent this issue is to reintroduce a length scale R over which the magnetic 
charge is smeared. This so-called “extended anyons” model is discussed in [MlEniEjj and is 
sometimes argued to be the correct physical description for anyons arising as quasi-particles 
in condensed-matter systems. In this paper we will allow R to become small when N —)• oo, 
in which case we recover the point-like anyons point of view, at least if one is willing to 
ignore the issue of non-commuting limits. 

Let us consider the 2D Coulomb potential generated by a unit charge smeared over a 
disc of radius R: 

wr{x) := log I . I * (x). (1.4) 

Observing that (with the convention wq := log | . |) 


V'’“r(;o(x) = -J— pr, with Bq{x) = V"*" • = A.wq = 2t:5q, 

|xp 

the natural regularization of Aj corresponding to the extended-anyons model is given by 

Af := ^ V^WR{xj - Xfc), (1.5) 


^For increased clarity we will in this work separate a from A, so that A corresponds to the statistical 
vector potential of fermions modeled as bosons. 
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leading to the regularized Hamiltonian 

N 

+ + V{xj)) . (1.6) 

i=i 

We shall denote 

E^{N) :=mfa{H§) (1.7) 

the associated ground state energy (lowest eigenvalue) for N extended anyons. 

For fixed ii > 0 this operator is self-adjoint on and one can even expand the 

squares to obtain a sum of terms that are all symmetric and relatively form-bounded with 
respect to the a = 0 non-interacting operatoi0. This gives 

N 

H§ = Y,ip]+yi^j)) 

j=i 

+ {pj ■ - Xk) + V^WR{xj - Xk) ■ Pj^ 

+ o? V^WR{xj - Xk) ■ V^WR{xj - Xi) 

+ a^'^\S/wR{xj-Xk)\^. (1.8) 

j¥^k 

We also note that by the diamagnetic inequality (see, e.g., [29l Theorem 7.21] for R > 0, 
and IMl Lemma 4] for R = 0) 

> (|T|,H^(a = 0)|T|), 

and hence E^{N) > NEq for arbitrary a, with Eq the ground state energy of the one-body 
operator Hi = p^ + V. 

1.2. Average field approximation. The few-anyons problem can be studied within per¬ 
turbation theory, yielding satisfactory information on the ground state and low-lying exci¬ 
tation spectrum gHuggiKio]. For many anyons however, it is hard to obtain results this 
way. A possible approximation to obtain a more tractable model when N is large consists 
in seeing the potential or dLSD as being independent of the precise positions Xj and 
instead generated by the mean distribution of the particles (whence the name, average field 
approximation [53]) 

A[p] := V-^wo * p, 

A^[p] :=V-^WR* p, (1.9) 

where p is the one-body density (normalized in L^(M^)) of a given bosonic wave function T 

p{x) = / \^{x,X 2 ,...,XN)?dX 2 ...dXN, 

Jr2(]V-1) 


By the boundedness of Vwr and using Cauchy-Schwarz, all terms are infinitesimally form-bounded in 
terms of Hn{ci = 0) and hence is a uniquely defined self-adjoint operator by the KLMN theorem |45l 
Theorem X.17]. We shall assume V is such that a form core is given by (7^(8^). 
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say the ground-state wave function. One then obtains from (11.61) the approximate A^-body 
Hamiltonian 

N 

^nIp] ■= {(p^ + NaA^lpiy + V{xj)'^ . 

i=i 

If one considers p as hxed, the ground state of this non-interacting magnetic Hamiltonian 
acting on is a pure Bose condensate 

'hv = 

where u G L^(M^) should minimize 

{u\{p + NaA%]f + V\u) = N-^ • 

For consistency, one should then impose that 

= P, 

which leads to a non-linear minimization problem. One thus looks for the minimum and 
minimizer of the following average-field energy functional (recall the notation fi ~ Na) 

N := (|(V + i/3A^[|u|2])u|%F|u|2) (1.10) 

under the unit mass constraint 



Note that, for this problem to be independent of N it is pretty natural to — in line with (|1.,3I) 
— assume that fi ^ Na is hxed. It is not difficult to see that if Na —>• 0 we recover at leading 
order a non-interacting theory, and we are back to the usual perturbation scheme. We should 
point out that the limiting functional which dehnes a strictly two-dimensional model 

of particles with a self-generated magnetic held B{x) = curl/3A[/?](x) = 27^(5p{x) without 
propagating degrees of freedom and to which one could further consider adding an external 
magnetic held, is also of relevance for various Chern-Simons formulations of anyonic theories 
(see, e.g., [SSI [Ml [13 E] ) • 

1.3. Average field versus mean field. In principle, the average held approximation does 
not require that the true ground state of be Bose-condensed. In fact, the most common 
application of it has been in perturbing around fermions a = 1 [131 [T9l ISHl ISTl [52] (this has 
even been argued to be preferable 131 ESI), and usually one even restricts to the homogeneous 
setting with p a constant. However, the case of hxed (5 ~ Na which is natural for the study 
of ()1.10l) . places the limit A ^ oo of the original many-body problem in a mean-held-like 
regime for bosons. Indeed, observe that in (|1.8p . the two-body terms in the second line 
and the three-body term in the third line weigh a total 0{N) in the energy in this regime, 
comparable to the one-body term in the hrst line. The two-body term in the fourth line is 
of much smaller order, 0(1) roughly, which is fortunate because of its singularity. Actually, 
if one takes bluntly R = 0, the potential |VrcoP appearing in this term is not locally 
integrable, and hence an ansatz Tat = would lead to an inhnite energy. For extended 
anyons, R > 0 and this term can be safely dropped for leading order considerations. 

The study of the regime (II.3p thus resembles a lot the usual mean-held limit for a large 
bosonic system (see [251 [261SZ] and references therein), but with important differences: 
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• The effective interaction is peculiar: it comprises a three-body term, and a two-body 
term which mixes position and momentum variables. 

• The limit problem (ll.lOli comprises an effective self-consistent magnetic field. A 
term in the form of a self-consistent electric field is more usual. 

• One should deal with the limit i? —)■ 0 at the same time as ^ oo, which is 
reminiscent of the NTS and GP limits for trapped Bose gases [521 IHTl IBUl 1^ HO] . 

In order to make the analogy more transparent, we rewrite, for any normalized A^-body 
bosonic wave function G 


iV-i = Tr 


N 


+ /3Tr Ipi- V-^wr{xi - X 2 ) + V-^wr{xi - X 2 ) ■ Pi] 


P) 




yN -2 

N-1 


Ti' 


(3) 


A^ - 1 


Tr 


V-^wr{xi - X2) ■ V-^wr{xi - X3 )) 
|VwR(a:i - X2)p7w^ 


( 1 . 11 ) 


where 

:= Tl'k+i^N [|^Ar)(^w|] 

is the /c-body density matrix of the state |'I'jv)(^Arl) normalized to have trace 1. The 
notation here means that we trace out the last N — k variables from the integral kernel of 

Since all terms at least at first sight weigh 0(1) or less, the folklore suggests to use an 
ansatz 


Ttv = 

= |u®^)(u®''l. (1.12) 

Inserting this in the energy, dropping the last term, which is of order N~^ at least for 
fixed R, we obtain to leading order 

N-^ {^n\H§\^n) ~ £r[u]. (1.13) 

Indeed, on the one hand. 


Tr 


(pi ■ V-^wr{xi - X2) + V^wr[xi - X2) ■ pi ) 

= i Vu{x)u{y) ■ V'^wr{x — y)u{x)u{y) dxdy 

— i u{x)u{y)'S/'^WR{x — y) ■ Vu{x)u{y) dxdy 

7 JR2x1R2 

= 2 / a^[|«| 2 ].JM, 


(1.14) 


using the definition (II.9p and denoting J[rt] the current 

J[u] := - {uVu — uVu). 


(1.15) 
















AVERAGE FIELD APPROXIMATION 


7 


Note that this is really a phase current density: 

J[m] = /oV(^ if u = 

On the other hand 


Tr 


V-^wr{xi - X2) ■ V-^wr{xi - X3)) 7^^ 


u{x)\'"‘\u{y)\'"‘\u{z)\'^'V'^WR{x — y) ■ V^wr{x — z) dxdydz 




(1.16) 


and it suffices to combine these identities in (jl.llh (and approximate [N — 2)/{N — 1) ~ 1) 
to obtain the desired expression ()1.13p for the energy. 


1.4. Main results. We may now state our main theorem, justifying the average held ap¬ 
proximation in the almost-bosonic limit at the level of the ground state. For technical 
reasons we assume that the one-body potential is conhning 

V{x)>c\x\^-C, s>0, (1.17) 

and that the size R of the extended anyons does not go to zero too fast in the limit iV —)• 00 . 
The rate we may handle depends on s. These assumptions are probably too restrictive from 
a physical point of view but our method of proof does not allow to relax them at present. 
Here and in the sequel, denotes the average-held functional (ll.lOp for 7? = 0, and 
its inhmum under a unit mass constraint. Although we do not state it explicitly, we could 
also keep R hxed when N ^ 00 and obtain the limit functional with hnite R. The case 
of anyons in a bounded domain is also covered by our approach (modulo the discussion of 
boundary conditions) and the results in this case can be obtained by formally setting s = 00 
in the following. 


Theorem 1.1 (Validity of the average field approximation). 

We consider N extended anyons of radius R ~ N~'^ in an external potential V satisfy¬ 
ing ()1.17p . We assume the relation 


0 < ?7 < r]o{s) 



(1.18) 


and that the statistics parameter scales as 


a = P/{N-1) 


for fixed /3 G M. Then, in the limit N ^ 00 we have for the ground-state energy 


E^{N) 

N 




(1.19) 


Moreover, if is a sequence of ground states for H^, with associated reduced density 
(k) 

matrices 7)^7 then modulo restricting to a subsequence we have 






( 1 . 20 ) 
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strongly in the trace-class when N oo, where yL is a Borel probability measure supported 
on the set of minimizers of , 

:= {u G L2(m2) ; 11^11^2 = 1, £^^[u] = 


Remark 1.2 (The size of the smearing radius). 

As we said before, taking R not too small in the limit —)• oo is a requirement in our 
method of proof. It is likely that localization arguments could allow to take an s-independent 
rj < r]o{oo) = 1/4, corresponding to the rate we obtain for anyons in a bounded domain. To 
obtain an even better rate would require important new ideas. 

It is in fact not clear whether some lower bound on i? is a necessary condition for the 
average field description to be correct. For very small or zero R, it is still conceivable that 
a description in terms of a functional of the form of (ll.lOp is correct in the limit. Indeed, 
our above restriction stems from the method used to bound the ground-state energy from 
below, while for an upper bound the conditions on R (and even the finiteness of /3) can be 
relaxed significantly. A further possibility would be to take short-range correlations into 
account via Jastrow factors, as in the GP limit for the usual Bose gas [321 EH EOl ESI SO]. 

It would in any case be desirable to be able to take R <C the typical interparticle 

distance in our setting, because one could then argue that smearing the magnetic charges 
has very little effect. Even in the formal case s = oo. Theorem 11.11 requires R S> S> 

which is rather stringent. We nevertheless obtain the functional for point-like anyons 
in the limit. 

It is sometimes argued in the literature [38l SOI E] that for anyons arising as quasi-particles 
in condensed matter physics, the magnetic charges should be smeared over some length scale 
R. The relevant relation between R and N then depends on the context. O 


The rest of the paper contains the proof of Theorem 11.11 We start by collecting in 
Section [2] some operator bounds on the different terms of the A^-body functional. This is 
required in order to have a correct control of the terms as a function of the kinetic energy 
in the limit i? —)• 0. For these estimates to be of use in the large-A^ limit we need an a 
priori bound on the kinetic energy of ground states of the A^-body problem, also derived in 
Section El We deal with the mean-field limit in Section El using the method of [26] . Some 
important adaptations are required to deal with the anyonic Hamiltonian, and we focus on 
these. The goal here is to justify (with quantitative error bounds) the sensibility of the 
ansatz when N becomes large, thus obtaining as an approximation of the 

ground state energy per particle. The basic properties of the average-field functional (ll.lOp 
are worked out in Appendix O In particular we study the limit i? —>• 0 to finally obtain 
as an approximation of the many-body ground state energy per particle. 

Acknowledgments. We thank Michele Correggi for discussions. Part of this work has 
been carried out during visits at the Institut Henri Poincare (Paris) and the Institut Mittag- 
Leffter (Stockholm). D.L. would also like to thank LPMMC Grenoble for kind hospitality. 
We acknowledge financial support from the French ANR (Project Mathosaq ANR-13-JS01- 
0005-01), as well as the grant KAW 2010.0063 from the Knut and Alice Wallenberg Foun¬ 
dation and the Swedish Research Council grant no. 2013-4734. 
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2. The extended anyon Hamiltonian 

In this section we give some bounds allowing to properly define and control the Hamil¬ 
tonian (HSl). As previously mentioned, for extended anyons, it is possible to expand the 
Hamiltonian as in (jl.8l) and estimate it term by term. By the boundedness of the interac¬ 
tion it follows that is defined uniquely as a self-adjoint operator on with the 

same form domain as the non-interacting bosonic Hamiltonian 

N 

+ V{xj)). 
i=i 

However, in order to eventually take the limit i? —>• 0 we will need to deduce more precise 
bounds depending on R. These will be used to deal with the mean-field limit in Section [3l 
In the following we introduce a fixed reference length scale Rq > 0, and always assume 
R Rq. Future constants, generically denoted by C, may implicitly depend on Rq. 


2.1. Operator bounds for the interaction terms. We start with some estimates on the 
different terms in (II.lip , exploiting the regularizing effect of taking R > 0. The following 
is standard: 


Lemma 2.1 (The smeared Coulomb potential). 

Let wr he defined as in dm). There is a eonstant C > 0 depending only on Rq sueh that 

sup IicrI < C-h I logi?|, sup IVwhI < sup iVtcnl < C*- (2.1) 

B(0,Ro) M2 K B{0,RoY 

Moreover, for any 2 < p < oo, 

( 2 . 2 ) 


Proof. A simple application of Newton’s theorem 
log |x| if \x\ > R 


wr{x) = 


log i? -I- 1 


|x|2 

'W 


— 1 I if 0 < |xl < i?. 


Theorem 9.7] yields 

x/jxp if jxj > R 
x/R^ if 0 < jxj < 


Vwr{x) = 


and ()2.ip clearly follows. For (12.2h we compute 

fR p fco 

II^'^«IIlp(m 2 ) = 271- + 2^ r-Prdr < R "P, 

where Cp > 0 depends only on p > 2. 


R, 

(2.3) 


□ 


We first estimate the most singular term of the Hamiltonian, corresponding to the fourth 
line of ()1.8p . Since it comes with a relative weight 0{N~^) in the total energy, the following 
bound will be enough to discard it from leading order considerations. 


Lemma 2.2 (Singular two-body term). 

We have that, as operators on T^(M'^) or 

\Vwr{x - 2/)P < C,R~^ {pI + l) 


for any e > 0. 


(2.4) 
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Proof. We start with a well-known simple application of Holder’s and Sobolev’s inequalities: 
for any W ; >-)• M and / G 

{f\W{x -y)\f) = [[ f{x,y)W{x -y)f{x,y)dxdy 

<\\W\\lp [ ([ dy 

JK2 VVR2 / 

<C\\W\\l, [[ {\VJix,y)\^ + \f{x,y)\^)dxdy 

J JR2xR2 

= C\\W\\L,{f\{-A, + l)®t\f) (2.5) 


where we may take any p > 1, q = G (1, -t-oo), and we use that in for any 1 < g < oo 

WgWl^q < Cq (^||V5||^2 -I- II 5 IIL 2 ) ) 

see, e.g., [29l Theorem 8.5 ii]. Next we may use (12.2p with W = iViUijp and p = 1 + e' to 
conclude 

WwWlp = W^wrWI,, < 

with a constant > 0 for given e > 0. □ 


We next deal with the two-body term mixing position and momentum, second line 
of (II.8p . This is somehow the most difficult term to handle, and it is crucial to observe 
that it acts on the current and not on the full momentum. We shall use three different 
bounds. In the following lemma, p2.6l) has a worst i?-dependence but it behaves better for 
large momenta than (12.7p and p2.8p . a fact that will be useful when projecting the problem 
onto finite dimensional spaces in the next section. Estimate (j2.8p might seem a bit better 
than (I22D, but we will actually need a bound on the absolute value in the sequel, which is 
not provided by ()2.8p . 


Lemma 2.3 (Mixed two-body term). 

For R < Rq small enough we have that, as operators on 

Px ■ V^wr{x -y) + ^^wr{x -y)-px < CR~^ \px\, 

Px ■ V-^wr{x -y) + V-^wr{x -y)-px < CeR~'' {pI -b 1), for all e > 0, 

and 


( 2 . 6 ) 

(2.7) 


± (px ■ V-^u;ij(x - y) + V-^WR{x - y) ■ px^ <(^(1-^1 logi7|) [pi -|- 1). (2.8) 

Proof. The bounds (j2.6p and (12.7p are based on the same basic computation. 

Proof of (12.6p . First note that 

Px ■ V-^wr{x -y) = V^wr{x -y)-px (2.9) 

because • V“'“r(;ij(x — y) = 0. We can then square the expression we want to estimate, 
obtaining 


[px ■ V^wr{x -y) + V^wr{x - y) ■ px'^ = ^Px ■ V-^wr{x - y)V-^WR{x - y) ■ px. 
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Consequently, for any / = f{x, y) G 

/| (px ■ V-^wjt(x -y) + V-^wji(x - y) ■ \f^ 


= 4 


Vxf(x, y) ■ V-^wr{x - y)) (Va;/(x, y) ■ V-^wr{x - y)) dxdy 


< 4 


Vx/(a:,y)|^ - y) dxdy. 


Inserting ()2.ip we get 


< 


_c 


K2xR2 


\Vxf\ dxdy 


and thus 


/ \2 (7 

■ V-^wr{x -y) + V-^wr{x - y) ■ p^^j < -^pI- 

We deduce (12.6p because the square root is operator monotone (see, e.g., [21 Chapter 5]). 

Proof of (j2.7p . We proceed in the same way but use Lemma 12.21 instead of just the rough 
bound (j2.ip (we denote x = (xi,a: 2 ) G 


/| (px • V-^Wi?(x - y) + V-^u;r(x - y )\f 


< 4 


// 

J Jr- 


|Va;/(®,y)|^ V-'-WRix-y) 


dxdy 


= ( dxj, 


V-^wr{x - y) 


dxj\ + '^-^WRix - y) dxj 

/ L2(r 4) \ / i2(R4) 

< ^ {{dxJ, {-Ax + 1) dxJ) + {dxJ, {-Ax + 1) dxJ)) < ^ (/, (-Aa, + 1)' /) 




Thus 


(^Px ■ V-^wr{x - y) + V-^wr{x - y) • px^ <^{pI + if 

for any e > 0, and the desired bound again follows by taking the square root. 

Proof of (12.8p . The idea is here a bit different. We pick / G C'“(ffi‘^;C) and compute as 

in (frni) 

(f\px ■ V^wr{x - y) + V^wr{x -y)-px\f') = 2 [[ V-^wr{x - y) ■ Ja;[/] dxdy 
' / 4 4r2xR2 

with 


ix[f\ = -{f'^xf-f'^xf). 





















12 


D. LUNDHOLM AND N. ROUGERIE 


We then split this according to a partition of unity y + 77 = 1 where x = 1 hr the ball 
B{0, Rq) and r] = 1 outside of the ball B{0, 2Rq): 

(f\Px ■ - y) + V-^WR{x - y) ■ Px\f \ = 2 // V-^ {x{x - y)wR{x - y))-3co[f]dxdy 


+ 2 


K2xR2 


V-^ {y{x - y)wR{x - y)) ■ J^if] dxdy. 


To control the x term we use Stokes’ formula and deduce 

2 // (x(x - y)wRix - y)yJM] dxdy = -2 // xix-y)wRix-y)cuTl,J4f] dxdy. 

It is easy to see that 

\cm%JM < 

pointwise, see e.g. HD Lemma 3.4]. We thus obtain 


± 2 

< 2 


ixix - y)wR{x - y)) ■ J^[f] dxdy 


IR2xR2 


\x{x - y)\\wR{x - y)\ \Vxf{x,y)\ dxdy 


< 


C{l + \logR\) [ [ \Va;f{x,y)y dxdy 

J 4R2xR2 

in view of (EH). For the y term we note that 

|Jx[/]|<|/||V,/|. 

Thus 

V-^ (y{x - y)wR{x - y)) ■ 3^[f]dxdy 


±2 


R2xR2 


< 


\^xf{x,y)\ dxdy 


V^{ri{x-y)wR{x-y)) \f{x,y)ydxdy + 

<R2 

<C [[ fi-A^ + 1 )/ dxdy. 

J jR'^xR'^ 

For the first term we used that (|2.3p implies 

y-^{ywR) < iVr/lIrcijj + jT/llVrciil < C 

because 77 = 0 in B{0,Ro) and 77 = 1 outside of i?(0,2iio)- Gathering these estimates we 
obtain the desired operator bound. □ 

The three-body term (third line of (ll.Sp l is actually a pretty regular potential term, as 
shown in the following: 

Lemma 2.4 (Three-body term). 

We have that, as operators on Lgyj^(M®), 

0 < V^wr{x - y) ■ V^wr{x - z)< C{pI + 1 ). ( 2 . 10 ) 

The essential ingredient of the proof is the following three-particle Hardy inequality of |18l 
Lemma 3.6] (see also |33] for relevant methods and generalizations): 
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Lemma 2.5 (Three-body Hardy inequality). 

Let d >2 and u : —)• C. Let TZ{x, y, z) be the eircumradius of the triangle with vertiees 

x,y,z & and p{x, y, z) := ■\/\x — y\^ + \y — z\^ + k ~ Then TZ~^ < 9p~^ pointwise, 
and 

3((i — 1)^ [ —}}^dxdydz < [ (\Vxu\‘^ + + \Vdxdydz. (2.11) 

jR3d p[x, y, zy J^3d V J 

Proof of Lemma Since we consider the operator as acting on symmetric wave functions 
it is equivalent to estimate 


Y] V-‘-wn(x - y) ■V-‘-wr(x - z). 

cyclic in x, y, z 

In general, let x,y,z € denote the vertices of a triangle and |x|r 
regularized distance. Then we claim the following geometric fact: 


0 < 


E 

cyclic in x, y, z 


x-y 

x-y\ji 


x — z ^ C 
X- zW ~ p{x,y,zy' 


( 2 . 12 ) 

max{|x|,i?} a 

(2.13) 


Recalling (j2.3p this gives a control on the expression we are interested in. Equivalently, we 
shall prove that 


0 < |y - z\ji{x - y) ■ {x - z) + \z - x\%{y - z) ■ {y - x) + \x - y\%{z - x) ■ {z - y) 

< (2.14) 

\x — y\ + \y ~ z\ + \z — x\^ 

for some constant C > 0 independent of R. 

Let us consider each of the different geometric configurations that may occur. In the 
case that all edge lengths of the triangle are greater than R, the cyclic expression that we 
wish to estimate in (|2.13l) reduces to ^TZ{x, y, z)~'^ (see [TSl Lemma 3.2]), which is clearly 
non-negative and bounded by |/0~^ from Lemma [2.5[ On the other hand, if all edge lengths 
are smaller than R then the expression equals -^^p^ (cf. [181 Lemma 3.4]), for which we 
have 0 < — ^P~‘^ since p^ < 3R^. If two of the edges are short and one long, say 

“ yli \y — z\ < R and \z — x\ > R, then the expression to be estimated in (|2.14l) reads 

R^{x — y) ■ {x — z) + \z — x\^ [y — z) ■ {y — x) + R? {z — x) ■ [z — y) 

' -V-' 

-{z-y)-{x-z) 

= R^ ((x -y) -(z- y))) ■ (x - z) + jx - zp {y - z) ■ {y - x) 

= \x — z\'^(^R^ + {y — z) ■ {y — x)^ > |x — z\^{FI? — R^) > 0. 

We furthermore have the upper bound 

\x — z\‘^ (^R^ + {y — z) ■ {y — x)^ < 2R^\x — z\^, 


while the r.h.s. of (|2.14l) is larger than 


R^|x-zp 
2R? + lx — zP 


> -R^\x — z\’^, 


using that |x — z| < |x — yj -|- ]i/ — 2 ] < 2R. 
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This leaves the case that only one edge is short, say \x — y\ < R, and the others long, 
\y — z\,\z — x\ > R. We thus consider the expression in (I2.14p 

\y - zp {x - y) ■ {x - z) + \z - xp {y - z) ■ {y - x) + R^{z - x) ■ {z - y). (2.15) 

We will here use methods from [33], namely the geometric (Clifford) algebra ^(M'^) over 
(see m for a general introduction). In the case d = 2 or d = 3 one can think of this as the 
real algebra generated by the Pauli matrices aj, with scalar projection (^)o := ^Tr^ and 
the embedding of scalars (0-vectors) 1 1 and of 1-vectors 9 x ^j'^j ^ 

and with the product of two 1-vectors xy = x ■ y + x A y decomposing into a traceful 
symmetric scalar part and a traceless antisymmetric bivector part. We have then, using 
tracelessness of the bivector parts of such products and the linearity and cyclicity of the 
trace, 

\y - zj^ (x - y) ■ (x - z) + jz - xp {y - z) ■ {y - x) 

= {{y - - y){x - z) + {z- xf{y - z){y - x))^ 

= {{y - z){x - y){x - z){y - z) + {y - x){z - x)(z - x)(y - z)}q 

= {{y - - y){z - x){z -y) + {x- y){z - x){z - x){z - y))^ 

= (((?/- -y) + {z-x){x-y) + 2(x -y) A{z- x)^ {z - x){z - y)^^ 

= {{y - x)ix - y){z - x){z - y))o 2 ((x - y) A (z - x) (z - x){z - y))^ 

= -{{x- yf{z - x){z - y))^ + 2{(x-^A{z - x) (z - x) A {z - y))^ 

z-y-{z-x) 

= -\x - 2/P {z-x)-{z-y) + 2{B^B)o, 

with i? := (z — x) A(z — 2 /) and its Hermite conjugate ill = (^z — y)A{z — x). In the fourth and 
fifth steps we used xy = yx + 2x Ay for the second term and then {y — z) + {z — x) = y — x, 
while for the final steps we again used the properties of the trace and that B^B = \B P is 
scalar. Thus, the expression (I2.15P we wish to estimate equals 

{R^ -\x - 2/P) (z - x) ■ (z -y) + 2|i?p > 0, 

where for the lower bound we also used that {z — x) ■ {z — y) >0 since x — y is the shortest 
edge. For an upper bound we can use permutation invariance (cf. [331 Proposition 15]) of 

|il| = |(x — 2/) A (x — 2;)| < |x — y\\x — z\, 

and for example that \y — z\ < i? -|- |x — z| < 2 |x — z\. Hence 

(i?^ — |x — 2/P) (z — x) ■ (z — y) + 2|ilp < 4 i?^|x — zp, 

while for the r.h.s. of (12.141) . with analogously \x — z\ <2\y — z\, 

R?\y — z\‘^\z — x\’^ R?\y — z\^\z — x\’^ 1 21 12 

R?+ \y — z\^+ \z — x\^ ~ 6|2/— zp 6 x z . 

We also remark that the non-negativity of (12.131) is in general false if | • jij is replaced by an 
arbitrary radial function, as can be checked when taking e.g. jxjij = 

Finally, the estimate (|2.1UI) follows simply by applying Lemma 12.51 with d = 2 to (I2.13P 
and using the symmetry of functions in Lgyj^(M®). □ 
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2.2. A priori bound for the ground state. For the estimates of the previous subsection 
to apply efficiently, we need an a priori bound on ground states (or approximate ground 
states) of the A-body Hamiltonian (11.611 . provided in the following; 


Proposition 2.6 (A priori bound for many-body ground states). 

Let S ® (sequence of) approximate ground states for that is, 

AT, at) < E^{N){1 + o(l)) when A —^ oo. 

Denote by 7 ^^ the assoeiated sequence of one-body density matrices. In the regime (11.311 . 
assuming a bound R > N~^ for some rj > 0 independent of A, we have 


Tr 


+ <C(1 + /32), 


(2.16) 


where C is a constant independent of (3, A and R. 


Proof. We proceed in two steps. 

Step 1. Using a trial state with u = \u\ ^ C'^(M^), we easily obtain from (ll.ll|l and 
the above bounds (note that the i?-divergent mixed two-body term is zero on such a u, and 
that the singular two-body term gives a lower-order contribution) 

E^{N) <C{1 +I3‘^)N. (2.17) 

Next we use the diamagnetic inequality [291 Theorem 7.21] in each variable to obtain 

N 

= V / (|(-iV,+aAf)T^|%U(x,)|Tjvp)dxi...dxAr 

^ 2M2iV V / 

N 

(|V,|Tjv||' + U(x,)|Tjvp) dxi... dxN. 


We deduce the bound 


Tr 


{p"+y)iNl <c{i+f3^), 


(2.18) 


where we denote 

7S,V:=Trfc+i^v[||'&iv|)(|'I^ivll] 

the reduced fc-body density matrix of ITat]. 

Step 2. Next we expand the Hamiltonian and use the Cauchy-Schwarz inequality for 
operators to obtain 

N 

^ (p^ apj ■ Af -h aA^ • pj a^\Af\‘^ -h V(xj)) 
i=i 

N 

> ((1 - 26-^)p] + (1 - 25)a^\Af\^ + V{xj)) 

i=i 


N 




R\2 
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choosing h = 4. Thus, using (I2.17D we have 


Tr 


+ <C(1 + /32) + 


C /32 


N 


N{N -ly 


\ i=i / 


(2.19) 


Then, since the last term in the right-hand side is purely a potential term 
IN \ / N \ 


Tjv, 

i i=i 


N ; = 


|T^l,J]|Af|2|T 
i=i 


AT 


We then expand the squares as in (jl.lip . and use Lemmas 12.21 and 12.41 to obtain for any 
e > 0 


1 


N 


N{N -ly l^f < CTr - xa) • V^n;R(xi - X3)7® 

+ CN-^ Th [|Vn;K(xi - 


< CTr 


(3) 


(pf -h 1) ® 1 (8) l7]v. 


-£ 7\r-l ' 


-h CeR~^N-^ Tr 


< C (1 + C,N-^R-^) Tr {pi + 1)7^;) 


( 2 ) 


(Pi + 1) ® ^7n. 


Inserting the estimate (|2.18p and recalling that we assume R> N ^ we conclude the proof 
by going back to (j2.19p . □ 


3. Mean-field limit 


We now turn to the study of the mean-held limit per se. The strategy is the same as 
in [26], but the peculiarities of the anyon Hamiltonian add some important twists, and we 
shall rely heavily on the estimates of the preceding section. 

3.1. Preliminaries. We hrst recall some constructions from |24l 126] . 

Energy cut-off. We denote by P the spectral projector of —A -|- V below 
energy cut-off A that we shall optimize over in the end: 

P := l/i<A, h = —A -|- V. 

Let 

Na = dim(PL2(M2)) 

be the number of energy levels obtained this way, and recall the following Cwikel-Lieb- 
Rozenblum type inequality, proved by well-known methods, as in |26l Lemma 3.3]: 

Lemma 3.1 (Number of energy levels below the cut-off). 

For A large enough we have 

Na < CA^+ 2/7 ^3 2) 


a given (large) 
(3.1) 


We shall also denote 


Q = l-P 


the orthogonal projector onto excited energy levels. 
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Localization in Fock space. We quickly recall the procedure of geometric localization, 
following the notation of [21]. Let 'Jn be an arbitrary Wbody (mixed) state. Associated 
with the given projector P, there is a localized state in the Fock space 

T{S)) = C © © • • • 


of the form 

off = G^ q © G^ i © • • • © G^ j^ © 0 © • • • (3-3) 

with the property that its reduced density matrices satisfy 




N 


n 


-1 N 

E 

k=n 


n 


Trn+l^.A: [G';A,fc] 


(3.4) 


for any 0 < n < N. Here we use the convention that 

7^^ := Tr„+i^Ar[7Ar], 


which differs from the convention of [2l| , whence the different numerical factors in (|3.4p . We 
also have a localized state G^ corresponding to the projector Q, which is defined similarly. 

The relations (13.4p determine the localized states G^,G^ uniquely and they ensure 
that G^ and G^ are (mixed) states on the projected Fock spaces F{PS)) and F{QS)), 
respectively: 


N 

E 

k=0 


Tr 


G 


P/Q 

N,k 


= 1 . 


(3.5) 


de Finetti measure for the projected state. We will apply the quantitative de Finetti 
Theorem in finite dimensional spaces of P El dZl [27] to the localized state G^, in or¬ 
der to approximate its three-body density matrix. The following is the equivalent of [261 
Lemma 3.4] and the proof is exactly similar: 

Lemma 3.2 (Quantitative quantum de Finetti for the localized state). 

Let'jM be an arbitrary N-body (mixed) state. Define 


dfiN{u) 


N 


E 



dfiN,kiu), dfiN,kiu) := dim{Pfi))')y^(^u'^^,Gff f,u‘^’‘'^du (3.6) 


and 


■= / |u®^)(n®^|d/iAr(ii)- 

JSPSj 

Then there is a constant G > 0 such that for every N G N and A > 0, we have 

GNa 


Tr 


pG3^(3)p®3 _ ~(3) 


< 


N 


(3.7) 


(3.8) 
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3.2. Truncated Hamiltonian. For an energy lower bound we are first going to roughly 
bound some terms in the Hamiltonian. Let us introduce the effective three-body Hamilton¬ 
ian 

'■=^{hi + h2 + hz)+ ^ ^ {pj ■V^WR{xj - Xk)+ '^^WR{xj - Xk) ■ 

^ l < jVfc <3 

+ I3^V^wr{xi-X2)-V^wr{xi-x^) (3.9) 

where hi is understood to act on the z-th variable (recall that h = — A -|- V). For shortness 
we denote 

W2=pi- V^wr{xi - X2) + V^wr{xi - X2) ■ Pi 
the two-body part of and 

VF3 = V^wr{xi - X2) ■ V^wr{xi - X3) 
its three-body part. With this notation 

H^-=l{hl+h2 + h3) + ^ W2{i,j)+/3^W3 

where W 2 {i,j) acts on variables i and j. Also note that for ||u|| = 1, by (|1.14p . (11.161) . 

=£l^[u] > Ef. 

We bound the full energy from below in terms of a projected version of H^\ 

Proposition 3.3 (Truncated three-body Hamiltonian). 

Let ^]\f be a (sequenee of) approximate ground state(s) for with assoeiated redueed 
(k) 

density matrices 'yffh Then, for any e > 0 and R small enough, 


— > Tr 


■^Rp®3^(3)pG3l +c'^ATr[Q7;;)i 


1 


a 


1 


^ iV + ^Rl+e + Ai?2 


Proof. We proceed in several steps. 
Step 1. We first claim that 


— {^N, Eft'll n) > Tr 


E3 In 


CgN-\ 


(3.10) 


(3.11) 


To see this, we start from (jl.lip . For a lower bound we drop the term on the fourth line, 
which is positive. Then one only has to correct the A-dependent factors in front of the 
third line. The term we have to drop to obtain (j3.11l) is bounded as 




N -2 
N -I 


Tr 


(v-^wr{xi - X 2 ) • \/^wr{xi - X 3 )^ 7 ^^ 


< C^N-^ 


upon using the a priori bound (12.161) combined with (I2.inp . 

Step 2. We next proceed to bound the right-hand side of (|3.11l) from below in terms of 
a localized version of and remainder terms to be estimated in the next step. We shall 
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need the projectors 

112 = 

IIs = 

and make a repeated use of the inequality 

ABC+ CBA>-eA\B\A-e-^C\B\C, e > 0, 


(3.12) 


for any self-adjoint operators A, B, C. 
We claim that 


Tr 


-^3 In 


> Tr 




Tr 


hQj‘}pQ 


|/3|(3 + 5i) Tt \u2\W2\Il2-f^i^] - |/3|<5r' Tr [p®2|^2|P®^J?^ 


(3) 


- 2|/3| Tr [P®" ® Q|1P2(1, 2)|P®" ® 

-(3\l + S 2 ) Tr [nslWglnsySJ^] - TT [p®3|lT3lP®37^3)] (3.13) 

where (ii and 62 are two positive parameters to be chosen later on. 

To prove (I3.13p . first note that 


P 3 In 

= Tr 

^7? 

+ - Tr 

2 


+ /3^Tt 



Tr 


Then, for the one-body term we have 


Tr 




= Tr 

> - Tr 
“ 3 


P/iPy^) 




+ Tr 

P®3 (hi + h2 + hs) P®='7j^"^ 


-F Tr 




using that h commutes with P and Q, PQ = QP = 0 and the fact that h is a positive 
operator. 

For the two-body term we write 


Tr 




= Tr 
-h Tr 


P®3lT2(l,2)P®37®j +Tr [n3lT2(l,2)n37{J^ 


jp®3w2(i, 2)n3 + n3W2(i, 2 )p®3) 7^) 

Next, since 

n3 = p®2 ® Q + n2 0 p + ns (g) Q, 
and 1 T 2 ( 1 , 2) only acts on the hrst two variables, this simplihes into 


(3.14) 


Tr 




= Tr 
-h Tr 
> Tr 


P®3iT 2(1,2)P®37®] +Tr [u5W2{1,2)U3jP 
{P®^W2{1, 2)U2 ® P + Hs ® P1F2(1, 2)P®3) 
'p^^W2{1,2)P^^jP] +Tr [n3lT2(l,2)n37S,'^ 


-5iTr 


(3) 


Hs (g) P|1T2(1, 2)|n2 (g) Pj}^ 


- <5f ^ Tr 


P®3|W2(1,2)|P®37® 
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where we use (I3.12p to obtain the lower bound. Then, using (I3.14p and (13.121) again for the 
second term of the right-hand side, as well as P, Q < 1 , we get 


Tr 


W2lf 


> Tr 




.,(2) 


- 2Tr 


-(3 + 5i)TV[n2|iT2|n27); 

Finally, the three-body term is dealt with similarly: 
Tr 


P^^ ^ Q\W2{1,2)\P^^ ^ Q-fP 


W^375?^] = Tr 

■(p®3 + n3)iF3(/^®^ + n3)7{?^ 



> Tr 


- (1 + 82 ) Tl 

n3|iF3|n37i,"^ 

-<^2-'Tr 

P®3|lT3|P®^7i^^ 


/oX 

using (I3.12P again. All in all, using also the symmetry of 7 ]^ , we obtain (I3.13|) . 
Step 3. We next estimate the remainder terms in (j3.13D . First we note that 


Tr 




> -TT 
“ 2 


Qi^n^Q 


+ — Tr 
20 




n37S^'^n3 


Va, 


+ —Tr 
4 


VhiU2j^^'^U: 


(3.15) 


The first inequality is just the definition of Q, and to see the second one we first write 


2 Tr 


VhQj^j^^Q = Tr VhiQ (g) -|- Tr 7 ^ 2 ! <8) Q 7 ® 


= Tr 
= Tr 
= Tr 
-I- Tr 
= Tr 
-h Tr 


-F Tr 


Vhi {QiSi P + Q^Q) 7^^^ 

Vhi{Q(E)P + Q0Q) + n2) 

\^2{p®Q + Q®Q) 7 ® {P®^ + n2) 
Vh2Q ® Qi^nQ ® Q 


Vh2 (P <Si Q + Q iSi Q) 7 ^^ 


^in27i,'^n2 


-|- Tr 


/i2 - '/hi) P (8) P ® Q 


> Tr 


VhiU2j^^'^n 


where we use repeatedly the cyclicity of the trace and the fact that i/h commutes with P 
and Q, along with the fact that PQ = QP = 0 and 

U2 = 1^^-P^^ = Q0Q + P0Q + Q0P. 

In the last step we also use that as operators 

'/h2P ® Q > a/AP ® Q > '/hiP (8) Q 
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by definition of the projectors P and Q. This gives the third term in the right-hand side 
of (I3.15p . The second one arises from similar considerations: 


Tr 




= Tr 

< 2Tr 
= 2Tr 

< 5Tr 

Next, using (I2.6p . we have 


n2 <8) Q'jP 


+ Tr 




+ Tr 


n2 <8) P'yP 


Q^n 

Qi^n 


+ Tr 
-h Tr 


" 27 }?’ 


{P(S)Q + Q®P + Q(S) Qh 


( 2 ) 

N 


Tr 


( 2 ) 


n2|Vh2|n27]v 




|pi|n27i'^n2 




Vhin27^^n2 


by operator monotonicity of the square-root, and by m 


Tr 


nsIfTsinsyi"^ 




n37j,'^n3 


(3.16) 


(3.17) 


Moreover, using 


again we get 


Tr 






P®^ ® QiP 


. cVa^ 


QIn'^ 


so that, combining with (13.1511 . choosing for some small fixed ci, C 2 >0 

( 5 i = ciVAR, 62 = C2AR? 

and A large enough (i.e. KR? > c for c large enough), we get 

Tr - l/3|(3 + <5i)Tr [n2|11^2|n27S,"^] - 2|/3| Tr ^ q|m/ 2 ( 1 , 2)|P®2 ^ q^(3) 

- /32(1 + 82 ) Tr [n3|iy3|n37}?^] > CTr [hQi^'^Q 
for some fixed constant C > 0. Then, inserting in (13.1311 . we deduce 


Tr 


^3 In 


> Tr 




TCTr 




C 


Var 

c 


Tr 


AP2 

(fc) 


Tr 


P®2|^^|P®2^(2) 

P®3|^^|P®3^J3) 


(3.18) 


But, using (|2.7p . ()2.16p . and Tryjy'^ = 1, 


Tr 


P®2|^^|P®2^(2) 


< CeR~^ Tr 


< CeR~^ Tr 


P®2(p2 ^ l)p®2^(2) 


( 1 ) 


{Pi + l)7iv 


<C,R- 
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whereas, using (12.10^ and (I2.16p again 


which completes the proof. 


< CTr 

< CTr 


P®3(p2 + 1)P®3^(3) 


[pI + l)7]v 


.( 1 ) 


<c, 


□ 


3.3. Energy bounds. In this subsection we prove the energy bounds establishing (I1.19p . 
The upper bound is obtained as usual by testing against a factorized trial state. Namely, 
taking in (|l.lll) with a normalized minimizer of and using (I1.14p . 

(jl.l 6 p , Lemmas 12.2112.3112.41 and the diamagnetic inequality (IA.3P , one finds 


K'l + (£r I“r1 +1) ( 


2 7J-£ 


N N 


= Pf + o(l) ^ (3.19) 


laf 


as R ^ N ^ with N —oo, where we also used Proposition IA.61 Note that for this upper 
bound we can allow any rate 0 < i] < oo, and may even take /3 = P{N) —)• oo. 

For the lower bound, inserting (|3.8I) in the estimate of Proposition 13.31 we get for any 
sequence of ground states 'if n that 

N 




+ CATr[g7}^)]-P^ 




-c 


1 


+ 


a 


+ 


1 


N VAPi+= 


> Tr 
-C 

We have here used Estimates 


"3 In 


CD A^+2/s 

+ CATr[Q^\i>] - (1 + I logPl) 


+ 


a 


+ 


1 


N VAPi+= 

and (I2.10p from Subsection 12.11 along with 
Pp^P < PhP < A 

to bound the operator norm of and Lemma l3.II to bound N\. 

Main term. Since by definition 7 ) 7 ' is a superposition of tensorized states we get 

Tr 


(3.20) 


P 3 In 


> Pf Tt 


~(3) 

Tn 


We then denote 


N 


A := TT 


= E 4 [Gr,* 


k=0 


(3.21) 


the fraction of P-localized particles. Using the simple estimate 


k k — 1 k — 2 k^ 
iVN- liV-2 ~ iw 


< CN-^ for 0 < A: < iV, 
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it follows from (j3.5p . (13.61) . (13.71) . and Jensen’s inequality that 


Tr 


~(3) 

7iv 


Since on the other hand 


we conclude 


I SPfi 


d^N = ^ 


^ 'N\-Wk 


k=3 


Tr 


N / u \ ^ 

^E(«) Tr [GJJ - 0(iV-') 

k=0 ^ ^ 

-I'l 


> A'^ - 0 {N-^) 


Tr 




= 1-A, 


Tr 


-^3 In 


(3.22) 


+ CATr[g 7 ^T] > ^3^af ^ _ , 7^-1 > ^af _ CN-^. (3.23) 


For the last inequality we bound from below in terms of the infimum with respect to 0 < 
A < 1. Since A is a very large number and is bounded as i? —)• 0 (see Proposition IA. 6 p . 
the infimum is clearly attained at A = 1. 


Optimizing the error. We next choose A to minimize the error in (I3.2np . We assume 
that R behaves at worst as 

N-'^. 


Changing a little bit rj if necessary we may ignore the | logi?| and R‘" factors, and we thus 
have to minimize 

1 1 A2+2/^ 

7lR + Ai^+ N ■ 


We pick 


A = _/vf^ 


to equate the first and the last term and get 


I I l^2{l+l/s) 

+Ai^+ N 


0(A^t(l+^)(l+V^)(l+4/5^)“^-l) + Q^jYtT+’?)(l+Vs)(l+4/5A"^-2)^ 


and this is small provided 


7 < % := 7 1 + - 


-1 


Since this is the main error term we conclude that the lower bound corresponding to p.l9l) 
holds provided R ~ with t] < rjQ, as stated in the theorem. The limit ^ E’^^ is 

dealt with in Appendix lAl 


3.4. Convergence of states. Given the previous constructions and energy estimates, the 
proof of p.2UI) follows almost exactly [2^ Section 4.3] and is thus only sketched. 

Modulo extraction of subsequences we have 
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weakly-* in the trace-class as N ^ oo. Prom Proposition 12.61 we know that (—A -|- 
is uniformly bounded in trace-class. Under our assumptions, (—A-|-U)“^ is compact and 
we may thus, modulo a further extraction, assume that 




( 1 ) 


strongly in trace-class norm. Then, by 


Corollary 2.4], we also have 
(fc) 


(fc) 

strongly for all A: > 0. 

Next we claim that the measure defined in Lemma 1,3. 2 1 converges (modulo extraction) 
to a limit probability measure /r G V{SS)) on the unit sphere of and that 


To see this, we first apply 
Tr 


7 ^^^ = / \u^^){u^^\dfj.{u) for all A; > 0. 

J uGSSj 

, with the above choice of A. We obtain 




ISPSj 




0 . 


(3.24) 


(3.25) 


On the other hand, combining (I3.23D with the energy upper bound (13.191) we get 

1 

where A is the fraction of P-localized particles defined in (13.211) . Using Jensen’s inequality 
as in (|3.22l) we deduce that 


fiNiSPSj) = TV rp®3^J^)p®3 
Combining with (|3.25p yields 


Tr 


(3) 

7iv 


ispr^ 


\u^^){u^^\dfiN{u) 


0 . 


(3.26) 


(3.27) 


Testing this with a sequence of finite rank orthogonal projectors 

K—^oo 

(3) 

and using the strong convergence of gives 

lim lim ijlm{Pk^) = 1, 

K^ooN^oo 

and we obtain the existence of a limit measure /r supported on the unit ball of by a 
tightness argument. Then (13.241) for k = 3 follows from (I3.27p . Since converges strongly, 
the limit has trace 1 and ji must be supported on the unit sphere. Obtaining (|3.24p for 
larger A: is a general argument based on (I3.26p . We refer to [26l Section 4.3] for details. 

There only remains to prove that /i is supported on But it follows easily from 

combining previously obtained energy bounds that 

dfi^iu) —)• 0. 

N^oo 


j 

Ei-ei[u] 

ISPS^ 
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Using in addition the results of Appendix]^ in particular Proposition IA.61 we obtain for a 
large but fixed constant C > 0 


/ 


S^f[u]<C 


_^afr 1 


djiNiu) —> 0 

N—>oo 


and 


/ dfiNiu) 0 . 

j£’^f[u\>C 

Then clearly /r must be supported on which concludes the proof. 


□ 


Appendix A. Properties of the average-field functional 


In this appendix we etablish some of the fundamental properties of the functional p.lUl) 
and its limit i? —)• 0. 

For S M and U : —)> M"*“ we define the average-field energy functional 


£ 


af 


[u] ■= (l (V -h ipA[\u\‘^]) u\^ + U|up) , 


(A.l) 


with the self-generated magnetic potential 

{x - y) 


A[p] := V-^tco *P = 


\x - y\ 


2 P{y)dy, curl A[p] = 27rp. 


The functional is certainly well-defined for u £ but we should ask what its natural 

domain is. We then have to make a meaning of for general u £ and the 

problem is that it is not certain that A[|up] £ even though u £ (see [12] for an 
exam pifl), so the product A[|up]tt may not be well-defined as a distribution (while Vtt 
certainly is). One way around this is to reconsider the form of the functional when acting 
on regular enough functions such that we can write u = \u\e^‘^ where (p is real. Then 

|(V -|- i/3A[|u|^])u|^ = |V|u| -|- i\u\{V(f + /3A[|rtp])|^ = |V|tt|p -|- ||rt|Vv9 -|- /3A[|rtl^]|ri||^ , 

where also Vp = \u\~'^^uVu and V|u| = |n| 
given by 

u 


-isp..-; 


3?uVu. Hence, an alternative definition is 


£^Uu\ := 


jR- 


iviuir + 


Vrt -|- /3A[|tt|^]|u| 


+ V\u\ 


(A.2) 


and the advantage of this formulation is that it makes clear that we actually demand 
|u| £ in order for < oo. We can then use the following lemma to see that in 

fact A[|ttp]rt £ L^(M^), and hence also Vrt £ L^(M^). (And conversely this also shows that 
if A[|ttp]ri ^ then we have no chance of making sense out of 

Lemma A.l (Bound on the magnetic term). 

We have for any u £ L^(M^) that 

[ |A[|up]|^ |up < ^||M||i2(ffi2) [ |V|u||^ . 


^Note that by Young’s inequality we have for any u G L^(R^) that A[|mP] G L^(R^) -|- eL°°(R^) for 
p G [1, 2). Also compare to the singular magnetic fields considered in |12II36| . 
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Proof. This follows from symmetry and from the three-body Hardy inequality of Lemma [T5 


A[|mP](x)|^ |tt(x)p dx = 

1 



x-y 
\x - 


X — z 


\x — z 


\u{x)\‘^\u{y)\^\u{z)\^ dxdydz 




n{xf 


\u\ 


^ix<\ 


|Vx| 


u 


<83 I 


'dX = 


If 


iVlrilr dx 


\u\‘^dx 


□ 


We can therefore define the domain of to be (and otherwise let := -|-oo) 


>af 


:= |u G : y < ooj , 

and we find using Cauchy-Schwarz, Lemma lA.il and |V|n|| < |Vri| that for u G 
0 < 8^^[u] < 2\\Vuf + 2p^\\A[\u\^]uf + j V\u\^ < (2 + 3P^\\uf)\\Xuf + J V\u\^ < oo. 
The ground-state energy of the average-field functional is then given by 


:=inf<iT“[d :u£ 


af r 


>af 


hr = 1 


For convenience we also make the assumption on V that V{x) —)• -|-oo as \x\ —)• oo and that 
C is a form core for ||n ||^2 := fjg ^2 with —A -|- V essentially self-adjoint 

and with purely discrete spectrum (see, e.g., [461 Theorem XIIL67]). This is then also a 
core for 


Proposition A.2 (Density of regular functions in the form domain). 

(^“(M^) is dense in w.r.t. namely for any u G there exists a sequence 

{Un)n^oo C C^(M^) such that 

0 o-'nd £^^[un] —^ as n ^ oo. 

Proof. Take u G , then ||Vtt||j ;^2 < oo and hence also ||u||lp < oo for any p G [2,oo) by 
Sobolev embedding. We use that C')T(M^) is dense in H^(R^), so there exists a sequence 
(^n)n—)-oo S.t. \\u-Un\\m ^0. Also, 

|||(V -7 f/3A[|rtp])u||2 - ||(V P il3A[\un\^])un\\2\ 

< II (V -I- il3A[\u\^])u - (V -b il3A[\unf‘])uh\2 

< ||V(n - Un )\\2 + I/3|||(A[|U|2] - A[\Un\^])u + A[\UnWu - Un)h 

< ||ri - UnWfji + |/3| ||A[|up - |MnP]tt ||2 + \/3\ ||A[|'UnP](tt “ U „)||2 , 
where by Holder’s and generalized Young’s inequalities 

llAd-up - |«„p]u ||2 < ||A[|up - |unp ]||4 Ihh ^ C \\\u\^ - \un?\\^/^ llVu^olla,^ II^^IU 
P C llzr — ^^lls/s — ^ 11^ — ^ ^5 

and similarly 

||A[|'ltjj| ] (rt 'lt,2)||2 Y a lilt ^nll/fl ^ 0) 

as n —)• oo. 

We also have continuity for ||n ||^2 here since we assumed that C'“(M^) is a form core. □ 
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Lemma A.3 (Basic magnetic inequalities). 

We have for u € that (diamagnetic inequality) 


|(V + il3A[\u\^])u^ > 


Vklp 


and 


|(V + i/3A[|ttp])'u| 



(A.3) 

(A.4) 


Proof. By density we can w.l.o.g. assume u € We then have A[|up] £ C 

-^foc(^^) hence the first inequality follows by the usual diamagnetic inequality (see e.g. 
Theorem 2.1.1 in mi)- Furthermore, by e.g. Lemma 1.4.1 in 


|(V + i/3A[|up])u|^ > ± 


curl (/3A[|up]) |up, 


which proves the second inequality since curlA[|up] = 27r|up. Instead of using density 
we could also have used the formulation (IA.2p or the fact that u G ^ A[|ttp] G 
p G (2, oo) by generalized Young. □ 


Proposition A.4 (Existence of minimizers). 

For any value o//3 G M there exists G with = 1 and . 

Proof. First note that for u G by Lemma lA.ll and Lemma I A. 31 

||Vtt||2 = ||Vtt + i/3A[|up]u — f/3A[|u|^]tt||2 < + |/3| ||A[|u|^]u||2 

< + W\^\\u\\l IIV|«|II2 < f 1 + l/ 3 |\/|l|.*ll^j 

Now take a minimizing sequence 

{Un)n^oo c ||u „||2 = 1, lim 

n—>-oo 

Then clearly (u„) is uniformly bounded in both L^(M^), Lfy, and (and hence in 

LP(M^), p G [2,oo)), and therefore by the Banach-Alaoglu theorem there exists G 
and a weakly convergent subsequence (still denoted Un) such that 

Un in l2(m 2) n n LP(Mf), Vun ^ in 

Moreover, since (—A + V + 1)“^/^ is compact we have that 

Un = (~A + F + 1) ^'^^( —A + F + V)^^'^Un 

is actually strongly convergent (again extracting a subsequence), hence 

Un u^^ in L^(M^). 

Also, A[|ttn|] converges pointwise a.e. to A[|mP] by weak convergence of Un in LP and, by 
the trick of Lemma lA.ll 
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by dominated convergence. The functions A[|unp]un are therefore even strongly converging 
to A[|ttp]u in by dominated convergence. It then follows that 

II(V + i/3A[|ttp])u|L = sup I(Vtt + i/3A[|up]tt, u)| 

lhll=i 

= sup lim \{Vun + iliA.[\un\^]un,v)\ 

ll^ll^in^oo 

< lim inf sup |(Vun + */3A[|unp]un, ^)| 

= lim inf ||(V + i/3A[|u„p])u„|| , 

n—>-oo " 


and since 11-11^2^ is also weakly lower semicontinuous (see, e.g., [SI Supplement to IV.5]), 
we have lim inf[u^] > Thus, with ||u®'^|| = limn->.oo ||rin|| = 1, we also have 


Proposition A.5 (Convergence to bosons). 

Let Eq resp. uq denote the ground-state eigenvalue resp. normalized eigenfunction of the 
non-magnetic Schrodinger operator Hi = —A + V, with V G L^^. We have 


E\ 




and that given an arbitrary sequence (tt,g) of minimizers for 


af 

(/3) 


UR —)• un in L 
^ /3^.0 


2/Tn)2\ 


up to a subsequence and a constant phase. 


Proof. Note that under our conditions for V, uq G is the unique minimizer of So = ^^(^^ 3 = 0 ) 

and can be taken positive (see, e.g., |29l Theorem 11.8]). By the diamagnetic inequality 
(]A.3I) . and by taking the trial state uq = |uo| in we find 

^0 < < q^)[uo] = So[uo] + /32 ||A[|uo|>o ||2 ^ 

(where we also used Lemma FA. Ih . and hence E^y.^ Eq as fd ^ 0. Now consider a sequence 

(u^) C of minimizers as /3 ^ 0 with —)• Eo, ||u,g|| = 1. Then, because of uniform 

boundedness and as in the proof of Proposition IA.41 we have after taking a subsequence 
that Up ^ u for some u G ||tt|| = 1, and also 

||Vu|| = sup KVu,u)| 
lhll=i 

= sup lim |(V'u,g+ i/3A[|u^p]u/3,u)| 

||t;|| = l P-^-O 

< lim inf ||Vu,g + z/3A[|tt,gp]u,g|| , 
so 

Eo < So[u] < £^^^[up]. 

It follows that So[u] = Eq and hence n = uq up to a constant phase. □ 
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From the bound (jA.4p we observe that the self-generated magnetic interaction is stronger 
than a contact interaction of strength 2-711/3| (despite the fact that we already removed a 
singular repulsive interaction in the initial regularization step for extended anyons). Hence 
we have not only > Eq by the diamagnetic inequality, but also 

E^^ > min f (27r|/3|/9^ -|- Hp) , (A.5) 

P>0, /b;2 P=1 Jm? 

which can be computed for given V by straightforward optimization. 

Let us now consider the corresponding situation for the regularized functional (extended 
anyons) 

M := (|(V + i/3A^[|np]) , A^[p] := * p, R>0. 

Since Vwr € L°°(M^) we have A'^[|up] € with 

iiA"[i<.nL < |hi 

and instead of Lemma lA.ll we have 

using Lemma 12.41 Hence the natural domain is again and all properties established 
above for are also found to be valid for (except ()A.4I1 and (IA.5D which now have 
regularized versions). Denoting 

E^ := min{£’^^[u] : u G ||tt ||2 = 1}, 

we furthermore have the following relationship: 


Proposition A.6 (Convergence to point-like anyons). 

The functional Ef^ converges pointwise to E^^ as R ^ 0. More precisely, for any u G 


Ei[u]-E^^\ 


u\ 


<C„|/3|(l + /34)(l + £:"f[u])3/2i?, 


(A.6) 


where Cu depends only on ||n|| 2 . Hence, 


E 


laf 

R 


R-j-0 


E 


af 


and if {ur)r^o C denotes a sequence of minimizers of Ef^, then there exists a subse¬ 
quence {uri)ri^o s.t. Ur/ —)• as R' —)• 0, where is some minimizer of E^^. 


Proof. We have for any u G that 

|||(V -t i/3A[|u|^])n||2 - ||(V i/3A^[|ttl^])u|y 

< II (V + i/3A[|u|2])n - (V + ipA^[\u\^])u\\^ = |/3| ||(A[|u|2] - A^dup])^^ 

< l/5| ||A[|uy - A^[|uy y ||n ||4 = 1/3| \\{Vwo - Vwr) * \u\^\\^ ||u ||4 , 
where by Young 

||(Vr;o - Vwij) * lupll^ < llVico - Vwijili |||wy|4 < ||Vu;o||ii(B(o,R)) hlls ^ 0> 
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as i? —)• 0, since Vwq G We deduce (IA.6D by combining this with previous 

estimates of this appendix and Sobolev embeddings. It follows that £fl[u] —>■ as 

R^O. 

Let {uji)R^o denote a sequence of minimizers of 

= II«rII = i, 

and take u G an arbitrary minimizer of . Then, since 

R-i-O 

we have that is uniformly bounded as i? —)• 0 and that 

limsupE'^^ < E^^. 

R^-o 

Then £^ [u/j], and hence also 

£^^\ur\ < C ^||ur||^i + II^rII^i^^ ^ C'^(^R [^r] +1)5 

are uniformly bounded as well. As in the proof of Proposition IA.4[ there then exists a 
strongly convergent subsequence {uri)ri^q, with ur/ —)• uq £ ■ Also, by weak lower 

semicontinuity T®'^[uo] ^ hminfjj/^o ^^^[rr']) so that for any e > 0 and sufficiently small 
R' > 0, 

^af ^ g^^luo] < £^^[uri] + e < £fl,[uR'] + 2e = E^, + 2e, 

where we also used that the convergence is uniform for our uniformly bounded sequence ur 
by the bound ()A.6D . It follows that < £^^[uq\ < E^^ + 3e, and hence uq is a minimizer 
with ||tto|| = 1 and E^^ = £^^[uo] = limR^oE^. □ 
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